Universal Approach to Optimal Photon Storage in Atomic Media 
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We present a universal physical picture for describing storage and retrieval of photon wave packets 
in a A- type atomic medium. This physical picture encompasses a variety of different approaches to 
pulse storage ranging from adiabatic reduction of the photon group velocity and pulse-propagation 
control via off-resonant Raman fields to photon-echo based techniques. Furthermore, we derive an 
optimal control strategy for storage and retrieval of a photon wave packet of any given shape. All 
these approaches, when optimized, yield identical maximum efficiencies, which only depend on the 
optical depth of the medium. 
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In quantum networks, states are easily transmitted 
by photons, but the photonic states need to be stored 
locally to process the information. Motivated by this 
and other ideas from quantum information science, tech- 
niques to facilitate controlled interactions between sin- 
glephotons and atoms are now being actively explored 
111, i, i, i, S, 0, 0, S, i, H El, E [H • A promising ap- 
proach to a matter-light quantum interface uses classical 
laser fields to manipulate pulses of light in opticall y d ense 
media such as atomic gases [2|, [sl, Ij, y, o, lD, Isl, Isl, llO|, luj 
or impurities embedded in a solid state material [iJ, ll3[ . 
The challenge is to map an incoming signal pulse into a 
long-lived atomic coherence (referred to as a spin wave), 
so that it can be later retrieved "on demand" with the 
highest possible efficiency. Using several different tech- 
niques, significant experimental progress towards this 
goal has been made recently [y, lZ|, |8| . A central question 
that emerges from these advances is which approach rep- 
resents the best possible strategy and how the maximum 
efficiency can be achieved. In this letter, we present a 
physical picture that unifies several different approaches 
to photon storage in A-type atomic media and yields the 
optimal control strategy. This picture is based on two 
key observations. First, we show that the retrieval effi- 
ciency of any given stored spin wave depends only on the 
optical depth d of the medium. Physically, this follows 
from the fact that the branching ratio between collec- 
tively enhanced emission into desired modes and spon- 
taneous decay (with a rate 27) depends only on d. The 
second observation is that the optimal storage process 
is the time reverse of retrieval (see also [ll', 'l35). This 
universal picture implies that the maximum efficiency is 
the same for all approaches considered and depends only 
on d. It can be attained by adjusting the control or the 
shape of the photon wave packet. 

A generic model for a quantum memory uses the A- 
type level configuration shown in Fig. [ija), in which a 
weak (quantum) signal field with freqeuncy v is detuned 



by A from the \g) — \e) transition. A copropagating (clas- 
sical) control beam with the same detuning A from the 
\s) — \e) transition is used to coherently manipulate the 
signal propagation and to facilitate the light-atom state 
mapping. In this system several different approaches to 
photon storage can be taken. In electromametically in- 
duced transparency (EIT) [2|, |3|, y, l5|, |6|, [D, ll2| , resonant 
fields (A = 0) are used to open a spectral transparency 
window, where the quantum field travels at a reduced 
group velocity, which is then adiabatically reduced to 
zero. In the Raman configuration [9, llO|, the fields have 
a large detuning (|A| ^ jd) and the photons are ab- 
sorbed into the stable ground state \s) by stimulated Ra- 
man transitions. Finally, in the photon-echo approach 
I 111 . ll3J . photon storage is achieved by applying a fast 
resonant 7r-pulse, which maps excitations from the un- 
stable excited state |e) into the stable ground state \s). 

A common problem in all of these techniques is that the 
pulse should be completely localized inside the medium 
at the time of the storage. For example, in the EIT 
configuration, a reduction in group velocity, which com- 
presses the pulse to fit inside the medium, is accompa- 
nied by narrowing of the transparency window, which 
increases spontaneous emission. Similarly, in the photon- 
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FIG. 1: (color online) (a) A-type medium coupled to a classi- 
cal field with Rabi frequency Q(t) and a quantum field with 
an effective coupling constant gyN . (b) Storage setup. The 
solid Q curve is the generic control shape for adiabatic stor- 
age; the dashed line indicates a vr-pulse control field for fast 
storage. For retrieval, the inverse operation is performed. 



echo technique, if a photon pulse is very short, its spectral 
width will be too large to be absorbed by the atoms. To 
achieve the maximum storage efficiency one thus has to 
make a compromise between the different sources of er- 
rors. Ideal performance is only achieved in the limit of 
infinite d [2]. 

In our model, illustrated in Fig. [ija), the incoming 
signal interacts with N atoms in the uniform medium of 
length L {z = to z = L) and cross-section area A. The 
control field is characterized by the slowly varying Rabi 
frequency Q{t — z/c). P{z^t) = V^Xli \9)i{^\/^zj where 
the sum is over all Nz atoms in a small region positioned 
at z, describes the slowly varying collective \g) — \e) co- 
herence. All atoms are initially pumped into level \g). As 
indicated in Fig.[T](b), we ffist map a quantum field mode 
with slowly varying envelope £in{t) (nonzero on t G [0, T] 
and incident in the forward direction at z = 0) to some 
slowly varying mode of the collective \s) — \g) coherence 
S{z^t) = v^^i\g)i{s\/Nz. Then starting at a time 
Tr > T, we perform the inverse operation to retrieve S 
back onto a field mode. As we explain below, the opti- 
mal efficiency is achieved by sending the retrieval control 
pulse in the backward direction; storage followed by for- 
ward retrieval is, however, also considered. The goal is to 
solve the optimal control [14] problem of finding the con- 
trol fields that will maximize the efficiency of storage fol- 
lowed by retrieval for given optical depth d = g'^NL/{jc) 
and input mode £in{t)- Here c is the speed of light, the 
atom-photon coupling g = p{v/{2heQAL)Y^'^ is assumed 
real for simplicity, and p is the dipole matrix element. 
The efficiency is defined as the ratio of the number of 
retrieved photons to the number of incident photons. 

Since the quantum memory operates in the linear 
regime, an analysis of the interaction process where all 
variables are treated as complex numbers is sufficient. In 
this limit the equations of motion read |l5| 



{dt^cdz)£{z,t) 

dtP{z,t) 



igVNP{z,t), (1) 

-(7 + iA)P{z, t) + ig^£{z, t) + 
MVt(t-z/c)S{z,t), (2) 

dtS{z,t) = in%t-z/c)P{z,t), (3) 

Here we have neglected the slow decay of S. For 
storage, the initial conditions are f(0,t) = fin(^), 
£{z,0) = 0, P{z,0) = 0, and S{z,0) = 0. Being 
the shape of a mode, £in{t) is normalized according to 
(c/L) Jq \£in{t)\'^dt = 1, so the storage efficiency is given 

by ^s = (1/^) Jo \S{z^T)\'^dz. For the reverse pro- 
cess, i.e. retrieval, the initial conditions are £{0^t) = 0, 
£{z, Tr) = 0, P{z, Tr) = 0, and S{z, T,) = S{L - z, T) for 
backward retrieval [16] or 5'(z, Tr) = S{z,T) for forward 
retrieval. The total efficiency in both cases is T^back/forw = 
(c/L) /^ \£out{t)\^dt, where £out{t) = £{L,t). 

It is instructive to ffist discuss the retrieval process. In 
a co-moving frame t' = t — z/c^ using a normalized co- 



ordinate C = z/L and a Laplace transformation in space 
C ^ 5, Eq. (dD gYve^£{s,t') = i^d-fL/cP{s,t')/ s. There- 
fore, the retrieval efficiency is given by 

il, = C-^Ud/{ss') j^ dt'P{sX)[P{s'\t')]*\ (4) 

where C~^ means that two inverse Laplace transforms 
(s ^ C ^^d s' ^ C,') are taken and are both evaluated at 
C = C' = 1- To calculate 77r, we insert £{s^t') found from 
Eq. P into Eq. ^ and use Eqs. ^^ to find 

dt {P{s,t') [P{s'\t't + ^(^,^0 [^(^",^0]*} 

= -7(2 + d/s + d/s')P{s, t') [P{s'\t't . (5) 

Eqs. (|4|5p allow us to express rj^ in terms of the initial 
and final values of the term inside the curly brackets in 
Eq. ([5j). Assuming P(s, oo) = 5'(s, oo) = (i.e. no exci- 
tations are left in the atoms) and taking £~^, we get 

7?r= / dC / dCka{C,C')S{C,T,)S*{C\T,), (6) 
Jo Jo 

fc.(C,C') = ^e-4i-K+^')/^)/o(Ml-C)(l-C')), (7) 

where /q is the zeroth-order modified Bessel function of 
the ffist kind. Note that rj^ does not depend on A and 
Vt{t). Physically, this means that a fixed branching ratio 
exists between the transfer of atomic excitations into the 
output mode ^out(^) and the decay into all other direc- 
tions. This ratio only depends on d and S{C,^T^). 

The efficiency t]^ in Eq. (|6]) is an expectation value 
of a real symmetric operator k^{CiC') in the state S{Q. 
It is, therefore, maximized when S{C) is the eigenvector 
(call it Sd{C)) with the largest eigenvalue r]^^^ of the real 
eigenvalue problem 



llrS{Q= / dC'fcd(C,C')5(C')- 



(8) 



To find Sd{C)i we start with a trial S{Q and iterate the 
integral in Eq. ([8|) several times. The resulting optimal 
spin wave Sd{^ — Q is plotted in the inset of Fig. [2] for 
d = 1, 10, 100, and d ^ oo. These shapes represent a 
compromise attaining the smoothest possible spin wave 
with the least amount of (backward) propagation. 

We now discuss storage. We claim that if, for a given 
d. A, and £in{t)^ we can find a control Q{t) that re- 
trieves backwards from Sd{l — C) into £^^{T — t), then 
the time reverse of this control, r^*(T — t), will give the 
optimal storage of £^in(^). To prove this, we represent 
our retrieval transformation as a unitary map U[Q{t)] 
in the Hilbert space H spanned by subspace A of spin- 
wave modes, subspace B of output field modes, as well as 
a subspace containing (empty) input and reservoir field 
modes (note that it is essential to include the reservoir 
modes, since the dynamics is unitary only in the full 




FIG. 2: (color online) Input mode ^in(t) (dashed) and control 
fields Q(t) (in units of ^JdrflT) that maximize for this ^in(t) 
the efficiency for resonant adiabatic storage (alone or followed 
by backward retrieval) at d — 1, 10, 100, and d ^ oo. Inset: 
Optimal modes 5*^(1 — C) to retrieve from backwards at d = 
1, 10, 100, and d ^ oo (C = ^/^)- These are also normalized 
spin waves S {(, T) / ^/rj^^^ in adiabatic and fast storage if it 
is optimized alone or followed by backward retrieval. 



Hilbert space of the problem). For a given unit vector \a) 
in A (a given spin wave), the retrieval efficiency is rj^ = 
\{b\U[n{t)]\a)\'^ = \{a\U-^[n{t)]\b)\'^, where we have used 
the unitarity of [/[l](t)], and where |6) is a normalized 
projection of U[^{t)]\a) on 5, i.e. the mode onto which 
the spin wave is retrieved. Introducing the time reversal 
operator T [H^e find r]r = \{a\rTU-^[n{t)]rT\b)\'^. 
One can show [l5| that the time reverse of the inverse 
propagator rU-'^[n{t)]T is simply U[n*{T - t)] so that 
we have r]r = \{a\rU[n*{T -t)]r\b)\'^. This immediately 
tells us that the time-reversed control Q^(T — t) will map 
the time reverse of the retrieved pulse into the complex 
conjugate of the original spin- wave mode with the same 
efficiency. The optimal spin waves are, however, real so 
that complex conjugation plays no role. Furthermore, the 
storage efficiency cannot exceed 77^^^^ since the time re- 
verse of such storage would then by the same argument 
give a retrieval efficiency higher than 77^^^^, which is a 
contradiction. Optimal storage is thus the time reverse 
of optimal backward retrieval and has the same efficiency 
^max _ ^max (^q^j^^ iuvolvcs the samc optimal spin wave). 
To identify the input modes, for which the optimal 
storage can be achieved, we use Eqs. ([T][3]) to analyti- 
cally solve the retrieval problem in two important limits: 
"adiabatic" and "fast". The "adiabatic" limit, whose 
two special cases are the Raman and the EIT regimes 
discussed above, corresponds to a smooth control field, 
such that P can be adiabatically eliminated in Eq. ([2]). 
Using the Laplace transform technique to eliminate £ 
from Eqs. (|l|2p . we reduce Eqs. (|2|3p to a simple differ- 
ential equation on S. We solve it, compute S, and take 
the inverse Laplace transform to obtain 



font Tr+-+t =-V^ dC^^e-^TT^ 

x/o(2y^dCMt)/(7 + iA))^(l-C,Tr), (9) 
where h{t) = J^ df\Q{f)\'^. We will now show that for 



a given d. A, and spin wave 5'(C), one can always find 
a control Q{t) that maps S{C) to any desired normalized 
output mode S2{t) of duration Tout, so that ^out(Tr+— + 
t) = ^/%£2{t) (provided we are in the "adiabatic" limit 
Toutdj > 1 [15]). To do this, we replace fout(Tr+^+t) 
in Eq. (|9|) with y^^2(^), integrate the norm squared of 
both sides from to t, change variables t -^ h{t)^ and get 



rj,J^dt'\£2it')f = ^CUhy^'dC 



7+iA 



-fdC + h' 
7 + iA 



xIo{2VWW/ij + iA))Sil-C,Tr)\, (10) 

which allows us to solve numerically for the unique h{t). 

Then |l](t)| = (^/i(t)) , while the phase is found by in- 
serting h(i) into Eq. ([9|). Optimal storage controls then 
follow from the time reversal argument above. Fig. [2] 
shows a particular Gaussian-like input mode E\^(i) and 
the corresponding optimal storage control shapes Vt pj[ 
for the case A = and d = 1, 10, 100, as well as the limit- 
ing shape of the optimal (^ as d ^ 00. As we have argued, 
the normalized atomic mode S{C,T)/^/rff^^^^ into which 
£in{t) is stored using these optimal control fields, is pre- 
cisely Sd{l — C)^ the optimal mode to retrieve backwards 
shown in the inset of Fig. [21 

The "fast" limit corresponds to a short and powerful 
resonant retrieval control satisfying Q ^ d'y that im- 
plements a perfect 7r-pulse between the optical and spin 
polarizations, P and 5*. This retrieval and the corre- 
sponding storage technique are similar to the photon- 
echo method of Ref. [ll|, ll3|]. Again using the Laplace 
transform technique, we find for a perfect 7r-pulse that 
enters the medium at time T^ 



jdL 



dCe-^^J^ 



/o(2y^)5(l-C,Tr), 



^out(Tr+-+t) = 

c 

(11) 

where Jo{x) = Io{ix). Since the fast retrieval control 
cannot be shaped, at each d, there is, thus, only one mode 
(of duration T ~ l/{jd)) that can be stored optimally. 
This mode is the time reverse of the output mode in Eq. 
(pT]) retrieved from the optimal spin wave Sd- 

We will now show that time reversal can not only be 
used to deduce optimal storage from optimal retrieval, 
but can also be used to find Sd in the first place. In the 
discussion above, the normalized projection of U~'^\b) on 
A (call it I a')) might have a component orthogonal to \a). 
In this case, the efficiency of U~^ as a map from B to 
A will be r][, = \{a^\U~^\b)\ > rjr- Now if the normalized 
projection of U\a') on B is not equal to |6), the map U 
acting on \a') will similarly have efficiency r]'J > ri[> rj^. 
Therefore, such iterative application of U and U~^ con- 
verges to the optimal input in A and the corresponding 
optimal output in B. Indeed, a detailed calculation J15j 
shows that the search for the optimal spin wave by iterat- 
ing Eq. ([8]) precisely corresponds to retrieving S{Q with 
a given control, time-reversing the output, and storing it 
with the time-reversed control profile. 




FIG. 3: ?7™ack (solid) and r]^^^ (dotted) are maximum total ef- 
ficiency for storage followed by backward or forward retrieval, 
respectively. ?7square (dashed) is the total efficiency for reso- 
nant storage of ^in(t) from Fig. [2] followed by backward re- 
trieval, where the storage control field is a naive square pulse. 



This time-reversal optimization procedure for finding 
the optimal \a) e A can be used to optimize not only 
retrieval, but also any map including storage followed 
by retrieval. For storage followed by backward retrieval, 
this procedure yeilds Sd{l — C) ^^^ maximum efficiency 
^bLck = ivT^'^)'^, since Sd{l — C) optimizes both stor- 
age and backward retrieval. Fig. [3] demonstrates that 
for resonant adiabatic storage of the field mode in Fig. 
[2] followed by backward retrieval, optimal controls result 
in a much higher efficiency f]^^^ than naive square con- 
trol pulses on [0,T] with power set by VgT = L (77square 
curve), where Vg = cVt^ /{g'^N) is the EIT group velocity 

i- 

For the case of storage followed by forward retrieval, 
iterations yield the maximum efficiency rf^^^ plotted in 
Fig. m It is less than rf^^^ since with backward retrieval, 
storage and retrieval are each separately optimal, while 
for forward retrieval a compromise has to be made. From 
a different perspective, forward retrieval makes it more 
difficult to minimize propagation since the excitation has 
to propagate through the entire medium. 

In conclusion, we have shown that the performance of 
EIT, Raman, and photon-echo approaches to a quantum 
light-matter interface can be understood and optimized 
within a universal physical picture based on time reversal 
and a fixed branching ratio between loss and the desired 
quantum state transfer. For a given optical depth (i, the 
optimal strategy yields a universal maximum efficiency 
and a universal optimal spin wave, thus, demonstrating 
a certain degree of equivalence between these three seem- 
ingly different approaches. We showed that the optimal 
storage can be achieved for any smooth input mode with 
Td^y ;» 1 and any A and for a class of resonant input 
modes satisfying Td'y ^ 1. The presented optimization 
of the storage and retrieval processes leads to a substan- 
tial increase in the memory efficiency. 

The results described here are of direct relevance to on- 
going experimental efforts, where optical depth d is lim- 
ited by experimental constraints such as density of cold 
atoms, imperfect optical pumping, or competing nonlin- 
ear effects. For example, in two recent experiments [6|,|7|, 
d ^ 5 was used, r]^^^ and r/square curves in Fig. [3] indi- 



cate that at this d, by properly shaping the control pulses, 
the efficiency can be increased by more than a factor of 2. 
Direct comparison to experiment, however, will require 
the inclusion of decoherence processes and other imper- 
fections. In Ref. [15], we discuss some of these imperfec- 
tions, as well as the details of the present analysis and its 
extensions to atomic ensembles enclosed in a cavity and 
to inhomogeneously broadened media. 

Finally, we note that the time-reversal based iterative 
optimization we suggest is not only a convenient mathe- 
matical tool but is also a readily accessible experimental 
technique for finding the optimal spin wave and optimal 
input-control pairs: one just has to measure the output 
mode and generate its time reverse. Indeed, our opti- 
mization procedure has been recently verified experimen- 
tally [l8|. We also expect this procedure to be applicable 
to the optimization of other linear quantum maps both 
within the field of light storage(e.g. light storage using 
tunable photonic crystals [l9[) and outside of it. 
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